function [t,T,T1]=main(T_set,v)   
    alpha=[0.004223340671227,0.003048958802801,0.001360914770792,0.002331976797780,0.001699887223936];
    h=[0.037615553555102,0.075031254321632,0.488031886585397,0.141783478826555,0.177666496009772];
    v=v/60;
    X=boundary();
    x=0:0.1:X(end);
    t=x/v;t=t';
    n=length(t);
    T0=zeros(size(t));
    for i=1:n
        T0(i)=get_K(T_set,x(i));
    end
        M=zeros(6,1);
    for i=1:4
        m=find(T0==T_set(i));
        M(i+1)=m(end);
    end
    M(6)=n;
    u0=25*ones(3,1);
    T=zeros(size(t));
    T1=T;
    for i=1:5
        [T(M(i)+1:M(i+1)),u0,u(:,M(i)+1:M(i+1))]=PDE(alpha(i),h(i),u0,T0(M(i)+1:M(i+1)),t(M(i)+1:M(i+1)));
    end
end
